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ABSTRACT 

Context. The magnetic topology of the solar photosphere in its quietest regions is hidden by the difficulties to disentangle magnetic 
flux through the resolution element from the field strength of unresolved structures. The observation of spectral lines with strong 
coupling with hyperfine structure, like the observed Mm line at 553.7 nm, allows such dilferentiation. 

Aims. To analyse the distribution of field strengths in the network and intranetwork of the solar photosphere through inversion of the 
Mn I line at 553.7 nm. 

Methods. An inversion code for the magnetic field using the Principal Component Analysis (PCA) has been developed. Statistical 
tests are run on the code to validate it. The code has to draw information from the small-amplitude spectral feature appearing in 
the core of the Stokes V profile of the observed line for field strengths below a certain threshold, coinciding with lower limit of the 
Paschen-Back elfect in the fine structure of the involved atomic levels. 

Results. The inversion of the observed profiles, using the circular polarization (V) and the intensity (I), shows the presence of magnetic 
fields strengths in a range from to 2 kG, with predominant weak strength values. Mixed regions with mean strength field values of 
1 130 and 435 Gauss are found associated with the network and intranetwork respectively. 

Conclusions. The Mn i line at 553 nm probes the field strength distribution in the quiet sun and shows the predominance of weak, 
hectoGauss fields in the intranetwork, and strong, kiloGauss fields in the network. It also shows that both network and intranetwork 
are to be understood at our present spatial resolutions as field distributions of which we hint the mean properties. 
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■ 1. Introduction 

' The absence of observables of the magnetism of the quiet sun is 
the main handicap in the determination of the magnetic topol- 
ogy of this region covering most of the solar photosphere. Over- 
, interpretation of the scarce information contained in the usu- 
. ally observed spectral lines in magnetometry, and its compari- 
' son with numerical simulations of magnetoconvection, has been 
the subject of long disputes about the nature of those fields. The 
addition of further observables, able to constrain the models and 
to unveil hidden ambiguities and biases in the diagnostic tech- 
niques was mandatory in such a situation. This work is a contri- 
bution in that direction, with a further step forward in the use of 
Mn lines with a strong hyperfine coupling. 

In a short descriptive summary of the previous studies of the 
magnetism in quiet Sun regions, we find from one side, the Fe 
lines in the visible region whose circular polarization amplitudes 
are almost uniquely sensitive to magnetic flux, (e.g. Keller et al. 
1994; Sanchez Almeida & Lites 2000; Lites 2002, Dominguez 
Cerdena et al. 2003; Orozco Suarez et al. 2007), and on the other 
side, we find the inversions of the Fe lines in the near-IR do- 
main whose profiles are almost universally split by Zeeman ef- 
fect (Lin 1995; Lin & Rimmele 1999; Khomenko et al. 2003). 
Traditionally, studies based on the visible spectral lines con- 
cluded on the presence of fields with kG strengths and small 
surface coverage, while studies based on infrared lines favoured 
weaker hG fields and more spread in surface. 
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Recently, the picture has been made more realistic by 
considering a continuum distribution function for magnetic 
field strengths at each resolution e lement in the quiet sun 
( Socas NavaiTO & Sanchez Almeida 2003) . instead of a single 
magnetic vector value. Additionally, a very interesting new ap- 
proach has enriched the individual line analysis in quiet regions: 
the inversions of simultaneous and co-spatial observations of 
the Fe lines in the visible (630 nm) and in the near-IR (1.5 
fim). Two such studies have already been caiTied out found- 
ing divergent results. In the first of these works, Dominguez 
Cerdena et al. (2006) have employed a model with three mag- 
netic components (mismas model, Sanchez Almeida & Landi 
Degl'Innocenti 1996) to simultaneously invert the lines in both 
spectral ranges. The authors retrieved field distributions in the 
second and third magnetic components with peaks in the strong 
regime and they conclude that the kG fields contribution dom- 
inates the magnetic flux and energy transport. In the other of 
these studies, Martinez Gonzalez et al. (2008) showed that when 
only the inversions of the Fe lines in the visible are considered, 
the magnetic distribution is dominated by strong strength fields 
(kG). However, when simultaneously near-IR and visible inver- 
sions were both considered, the kG contributions disappears and 
the distributions are mostly dominated by weak strength fields 
(hG). This last conclusion is in accord to the results obtained 
from near-IR inversions alone (see previous references). These 
results, suggesting a probable bias of the quiet sun Fe inver- 
sions when performed only in the visible at 630 nm had also 
been found on numerical tests (iBellot Rubio & Collados 20031 ; 
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iMartmez Gonzalez et al. 2006h . In the counterpart, a possible 
bias in the inversions of the IR alone, has been also argued sug- 
gesting that this lines are slightly sensitive to the strong strength 
fields, such that when mixed and not resolved magnetic struc- 
tures of hG and kG are both present, the inferred strengths of 
hG are preferred over those of kG (e.g. Sanchez Almeida & 
Lites 2000; Domfnguez Cerdena et al. 2006). Such discussions 
and differing results illustrate the difficulty and ambiguity of the 
measurements attempted, and have prompted the search for more 
observables. 

Apart from those Zeeman-based techniques discussed above, 
the description should be completed with the Hanle effect diag- 
nostics, pioneered by Stenflo (1982). A paradigm of not struc- 
tured and turbulent fields, e.g. (Manso Sainz et al. 2004), could 
result in apparent contradiction with the tube-like structures used 
to describe Zeeman-based observations if one insists in the pic- 
ture of a single vector magnetic field per resolution element. 

Mn lines subject to strong hyperfine coupling appeared as 
an interesting observable to be added to the previous ones: their 
amplitudes are dependent on magnetic flux as the usual Fe lines 
in the visible regime, but under certain magnetic field regimes a 
new spectral feature appears in the Mn profile. Just the presence 
or absence of such a spectral feature allows the observer to de- 
termine the presence or absence of fields with strengths below a 
certain threshold thus adding a new observable on field strength 
to be compared to measurements with near-IR lines. 

Quiet sun conditions limit the information on magnetic fields 
to circular polarization profiles, whose amplitudes stay too near 
to the usual signal-to-noise levels found in solar polarimetry. The 
amplitudes of the linear polarization Stokes parameters use to 
be below the noise level at the typical spatial resolutions, and 
are thus useless for analysis purposes; and intensity profiles are 
insensitive to the weak magnetic fluxes involved. Consequently 
one does not only need to add further spectral lines to the amount 
of observables, but also do it under observing conditions that 
guarantee very high signal-to-noise levels in polarimetry. We 
used the telescope THEMIS for that purpose in this work and 
we explored the capabilities of the Mn lines to provide further 
information on the quiet sun magnetic strength distribution. We 
simplified the formation of the Mn lines with a radiative transfer 
model using a Milne-Eddington atmosphere; the model proves 
to be at the limit of present possibilities with yet too many of its 
parameters left undetermined by the available observables. We 
nevertheless prove through synthetic tests that we can ascertain 
the distribution of field strengths in strong and weak flux regions 
(photospheric network and intranetwork respectively) with good 
accuracy, within the limits of the adopted Milne-Eddington mod- 
eling for the Sun. 

With the results presented here, that show a field strength 
distribution with maximum number of occurrences at very weak 
fields (less than or around 100 G), we propose in the conclu- 
sions a cartoon scenario that would be consistent with both 
the Zeeman-based measurements (visible and near-IR) and the 
Hanle measurements. 

In section 2, we review the signatures of the magnetic field 
strength in the line formation of the Mm profiles. In section 3, 
we describe the inversion code used in this work to retrieve the 
strength and magnetic flux from the data, and we test its ca- 
pabilities with realistic noise conditions. The observations are 
described in section 4 and the results presented in section 5. 
Conclusions and a final overview are included in the final sec- 
tion. 



2. Hyperfine structure and the Stokes parameters 

The coupling of the nuclear angular momentum (/) with the or- 
bital total angular momentum (J^s ) gives rise to hyperfine struc- 
ture (hfs). The total angular momentum F - J + I results in 
new non-degenerated atomic levels. The apparent splitting of the 
atomic transitions due to this effect, and the consequent broad- 
ening of solar and stellar spectral lines has been discussed in 
several papers since the work of Abt (1952), and used to mea- 
sure the abundances of different elements (Kurucz 1993). 

The most common observable for the analysis of (solar 
or stellar) spectral lines modified by hyperfine structure is the 
broadening of the line profile. Such broadening however can 
be the result of others physical phenomena as well: micro- 
turbulence, Doppler broadening or just blended lines. Too of- 
ten the ambiguity on the origin of an observed line broaden- 
ing cannot be solved for, even in observations with very high 
spectral resolution, unless full spectropolarimetry is performed. 
Such is the case presented in this work, the Mm line at 553.7 
nm: although producing a slightly asymmetric intensity profile, 
the recognition of the hfs features becomes straightforward in 
circular polarisation, since the hps together with an atmosphere 
permeated of an external magnetic field produce a unique spec- 
tral signature in the Stokes V profiles. The use of this signa- 
ture and its variations under different conditions as a tool for 
the detection of magnetic fields in the quiet sun in different field 
strength regimes has been theoretically first proposed in Lopez 
Ariste et. al (2002). Later, in Lopez Ariste et. al (2006), the au- 
thors showed that observationally it is possible to detect the hfs 
effect in the circular polarization profiles. In this work we extend 
previous studies and we build an inversion code drawing its in- 
formation from those spectral features in the Stokes profiles due 
to HPS. The inversion algorithm is firstly validated and charac- 
terized with synthetic data, and then applied to observed data. A 
quantitative analysis of the field strength distribution in the quiet 
sun results. 

For illustration and reference we show, in Fig.[Tl the calcu- 
lated spectral signature of the Stokes V profiles for the Mn i at 
553.7 nm, as a function of the magnetic field strength. The im- 
portant remark from the V profiles pictured is that the weaker the 
magnetic field strength, the higher the amplitude of the spectral 
feature due to the hps, visible in the central part of the profile. 

3. The inversion code 

3.1. Solution to the radiative transfer equations 

The four coupled equations for the radiative transfer of the 
Stokes parameters were solved analytically by the code diago- 
nal (Lopez Ariste & Semel 1999), which was slightly modified 
from its original version to include the computation of hyperfine 
structure. 

To use the code diagonal we assume a Milne-Eddington 
atmosphere, with the possibility to include numerically gradi- 
ents with optical depth of the velocity and magnetic fields, and 
adding up to fifteen free parameters. Evidently, the use of a big- 
ger number of free parameters implies a higher risk of inversion 
ambiguities and biases, if they happen to exceed the dimension- 
ality of the observables. But on the other hand considering less 
parameters may not allow to cover all the variety of the observed 
profiles. We decided thus to limit the number of parameters for 
each atmospheric model, to the ones described below, leaving 
out the velocity and magnetic gradients along the line of sigth. 

The inversion code works in hierarchical mode by a double 
pass fit. In the first pass only the intensity is fitted, permitting 
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Fig. 1. Stokes V profiles as a function of the magnetic field 
strength. Note that the spectral feature on the central part of the 
profiles is the signature associated to the hps regime, and that 
it disappears about B ~ 750 Gauss, as the Une passes into the 
incomplete Paschen-Back regime. 

half of the atmospheric parameters (Doppler width, Yo, velocity 
in the line of sight Ylos^ and the source function gradient, 
Vt) to be constrained in their ranges of variation. The remain- 
ing four atmospheric parameters are related to the presence of 
a magnetic field: the angle between the magnetic field direction 
and the direction of the line of sigth Oblos , the azimuthal angle 
0;, the strength of the magnetic field B and a filling factor //, 
accounting for unresolved structure. In the second (and final) fit, 
the eight parameters are considered in fitting simultaneously the 
profiles of intensity (I) and circular polarization (V). Each atmo- 
spheric parameter will be left to vary randomly in the range of 
values specified in table [T| in order to build a database for the 
inversions. 



3.2. Atmospheric models and database 

Let X denote a given combination of the described atmospheric 
parameters such that. 



Xm = {Vd, Vlos , ri(h Vt, B, Oblos , OzJf),, 



(1) 



is the m,k combination. Let S„,(Xm) denote the Stokes vector cal- 
culated for this combination of atmospheric parameters 



(2) 



and let P,„ denote any of the components of S„,{Xm). 

In the computation Stokes parameters we have considered 
201 spectral wavelength points with a step size of 12 niA, cover- 
ing in total 2.41 A around the line center at 553.762 nm. In the 
following, we proceed to consider each profile P,„ as a vector of 
201 points P„,. The creation of the database is the most impor- 
tant step for the inversions working with principal components 
analysis (e.g. Rees et al. 2000; Socas-Navarro et al. 2001). It is 
required to consider a large enough' number of profiles in or- 
der to consider the database as statistically representative. Thus, 
firstly, we have created a database of 80804^ profiles P,„, here- 
after called the profile space matrix (Ai), and we assume that all 



' The actual figure has no particular meaning other than illustrating 
what we consider as large enough and for future reference. 



possible observed profiles in the data set are represented in the 
database. This assumption is however verified in the tests pre- 
sented at the end of this section. 

We next perform a S ingular Value Decomposition (SVD) 
jGolub & Van Loan 19961) of the space matrix to compute an al- 
gebraic basis of the space spanned by the profiles. We should 
refer to the elements of the basis as eigenprofiles. 



3.3. Analysis with PC A 

Following the synthetic notation of Eq. for each one of the 
Stokes parameters (/, Q, U, V) it is possible to construct the 201 x 
80802 profile matrix M that contains all the profiles P,„ of the 
database: 



M = (Po 



80801 



(3) 



After computation of the eigenprofiles, each profile in the matrix 
Ai can be reconstructed as: 



PmiXm) = ai(jc„) Ci ; / = 0, 1, 200 



(4) 



where c, are the eigenvectors associated to the space profile {Ai) 
and a, are the coefficients of the linear combination associated to 
the profile P„,. Since the eigenvectors form an orthonormal basis 
of the space (Al) arranged in order of importance, one can cut 
the expansion in Eq. (HJi and preserve only the first eigenvectors 
that contain relevant information on the shape of the profile. The 
rest of eigenvectors are asso ciated with the noise and there i s no 
interest in preserving them jHansen 19921: iRees et al. 20001) . In 
our study we find that the appropriate number of components in 
Eq. (Ell, is i,„ax - 30. 

With the help of the set of eigenvectors, {e,), we create a 
second database {Ai') that is independent from the {Ai) matrix, 
used to compute the {c,}. The {Ai') database will be decomposed 
in terms of the coefficients a, as follows. Let us denote P„ a 
vector of the new (second) database. Using the Eq. (|4|, 



(yi(Xn) = Pn(Xn) ■ ei ; / = 0, 1 , 29. 



(5) 



In this way, each profile P„ (observed or calculated) has asso- 
ciated a unique set of coefficients (ao, ... , awmx), denoted here- 
after a„ = {oq, ... ,ai„,ax). So instead of working in the profile 
matrix space (Al') it is equivalent, and faster, to work in the co- 
efficient's space (C). From Eq. (O, and in analogy with (Al'), 
we obtain for each profile P„ the correspondent coefficients or,, 
in order to build the 30 x 80802 coefficient matrix (C), 



C = (ffo • ■ ■ am%o\f ■ 



3.4. Description and validation of the inversion code 



(6) 



One important quality of the code developed is that it is very 
easy to estimate the error bars per model parameter. We present 
here the tests applied for that purpose before pursuing into the 
inversion of the observed profiles. 

Let Psy„ denote a synthetic profile with associated coeffi- 
cients cCj,,,. The profile considered to be the solution by 
the inversion code, is found in the space of coefficients (C) as 
the nearest set of coefficients (in an Euclidean sense) to asy„ : 



\\ai,n.-a,„.„\\ = min{\\asy„-a„\\) ; « = 0, 1, 80801. 



(7) 



Since each one of the Stokes parameters would give, if 
considered separately, one solution to the Eq. (|7]i, let us say 
(Qr,„,,,/, a,>,i,,g, dinwUi cuinv.v), it is required to find the best solution 
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Table 1. All the parameters vary randomly in the indicated ranges. The first four parameters of the initial set, are used to fit the 
intensity in the first inversion passage in the code. In the final set, all the eigth parameters are considered to invert simultaneously 
(I) and (V). 



Parameters 


Vd (mA) 


10 


Vr 


Vlos (km/s) 


S(G) 


Oblos (deg.) 


0z (deg.) 


ff 


Initial Set 
Inverted (max,min) 
mean value 


[35,75] 
(41,60) 
47 


[0,2.5] 
(0.7,1.8) 
1.16 


[0,15] 
(1.0,14.9) 
8.8 


[-4,4] 
(0,0.54) 
0.02 


[0,2000] 


[0,90] 


[0,90] 


[0,1] 


Final Set 


[39,61] 


[0.5,2] 


[0,15] 


[-2,2] 


[0,2000] 


[0,90] 


[0,90] 


[0,1] 



for the four Stokes parameters simultaneously. For this purpose 
and given that it is a well-known fact, but particularly true in 
our data sets, that /„ » V„ > U„, Q„ necessarily implying that, 
in average, a,-,,,,,/ > Qrm,.,v > cuim-.u^ ctim;Q (e-g- Lopez Ariste & 
Casini 2002), we find useful to introduce a weight for each solu- 
tion. Thanks to this, the polarized signals are wighted similarly 
to the intensity one, despite their disparate amplitudes. 
Let the scalars 



max(a synj) 
max{asyn,s ) 



■,S ^I,Q,U,V 



(8) 



be the factors that modify the computation of the norms in Eq. 
(|7]), such that each one of the four solutions has the appropriate 
weight and can be added to the others in a linear way. The final 
solution, asoi, is then found as: 

{l.s aiv.s I \amv,s - 0!sy„,s 1 1) = min {l.s a^,,s I \a!sy„,s - 1 1) (9) 

where S denotes, as usual, the Stokes parameters (I,Q,U,V). 

Since we know from observations, described in section 4, 
that the signal-to-noise ratio in linear polarization (Stokes Q, U) 
is not big enough to include them in the analysis, they were ex- 
cluded as well from the tests applied to the code and only the 
intensity and circular polarization were considered. This implies 
that from Eq. (O, 5 = /, V. 

Finally, we compare the solution values of Xsoi (associated 
to asoi) to the known values of Xsvn (associated to asy„). This 
comparison leads to the statistical determination of the inversion 
errors per model parameter, where by error we mean the absolute 
difference between the original and the retrieved values. 



3.4.1 . Discussion of tine tests 

In order to not bias the tests we have constrained the thermody- 
namic parameters to the ranges of variation retrieved from the 
observational data sets after the first inversion (see the first three 
rows of table [T]i. Otherwise, we have not imposed any a priori 
constrain to the magnetic parameters, and the profiles used in the 
tests were synthesized using independent uniform distributions 
for each parameter, following the values found in the last row 
of table [T] However, in the search of direct translation to the in- 
version of quiet sun profiles, we have conditioned the profiles to 
have amplitudes in the V profiles that do not exceed the maxi- 
mum amplitude found in the data set (V < 5 x 10"^) while being 
bigger than the lower noise level (V > 5 x 10"^). Under these 
conditions, it is expected that the sample of profiles used in the 
tests is comparable with that from quiet sun data. 

In Fig.|2l we show the results of the inversion of samples of 
500 profiles. As a first test, we have not included noise in the 
inverted profiles. With the exceptions of Vr and 6^, the errors in 
most of the parameters decrease rapidly to zero. 

The bad results for 0, are due to the fact that the linear polar- 
ization. Stokes Q, U, was not included in the inversions, imply- 
ing that the inverted values in this parameter correspond to an 
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Fig. 2. Inversions histograms of each atmospheric parameter. 
The number in the upper right corner indicates the value of the 
bin size used. 



aleatory distribution. Consequently, the original (by construc- 
tion) and the inverted distributions in the azimuthal parameter 
are aleatory. Finally, plotting the absolute difference between 
these two distributions results in an histogram with a pick er- 
ror centered in (around) zero and that decrease smoothly, like 
the one obtained for 0,. 

Now, the most relevant parameter in our analysis is the mag- 
netic field strength B, since the major interest of this work is to 
analyse the magnetic field in quiet sun conditions. In what fol- 
lows we focus the discussion on this parameter and we include 
also the results of inversion of profiles with noise. 
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3.4.2. Analysis of effective inversions in B 

The graphics in Fig. |2] correspond to the ideal case of profiles 
without noise. Considering that the real profiles have a typical 
noise level of the order of 2 x 10 to 5x10"^ the continuum level, 
we proceed to invert synthetic profiles with noise added to those 
levels. Let B,„i, denote the retrieved value of the field strength, 
and let B ,,,, denote the real atmospheric value of the synthetic 
profile. Representing by (AB) the upper limit of the error value, 
such that 



(10) 



we can characterize the inversion efficiency behaviour of the 
code, understood as the percentage of profiles correctly inverted 
as a function of (AB). In tabled we compare the inversions of 
the magnetic strength for profiles with different noise levels, and 
in Fig.[3]we plot the respective cumulative distribution functions. 



Table 2. Inversions efficiency of the magnetic field strength, ex- 
pressed in percentage of correct inversions as function of AB . 



AB 




Noise level 









5x10-5 


2x10-* 


50 G 


69 % 


46 % 


40 % 


100 G 


89% 


74% 


65 % 


150 G 


96 % 


88% 


80% 


200 G 


99% 


93 % 


89% 


300 G 


100% 


98% 


96% 



c 
o 



CD 
> 



o 

o 
o 




No Noise, AB=150G, 96% 
Noise = 5xl0"', AB=150G, 88% 
Noise = 2xl0"*, AB=150G, 80% 



100 200 
A B (Gauss) 



300 



Fig. 3. Curves of cumulative percentage of correct inversions. 
Each sample consisted of 500 profiles. The solid line corre- 
sponds to profiles without noise. The dashed line corresponds 
to a noise level of 5x10-^ and the dash-point-dash line to a noise 
level of 2 X 10-"* the intensity of the continuum. The inset indi- 
cates the efficiency percentage for an upper limit of AB =150 
G. 



The value of the magnetic field strength goes from through 
2000 G covering all the range of magnetic strengths inherent 
to the quiet Sun regime. Since the spectral signature in Stokes 
V profiles evidence of hps regime, like those in Fig. [1] disap- 
pears for values around B ~ 750 G, we have special inter- 
est in the behaviour of the code while inverting profiles in the 



weak field regime. In table [3] we show that the precision invert- 
ing the magnetic strength is mainly independent of the magni- 
tude of the magnetic field. The inversion gently improves with 
higher values of the field strength, but in general the precision 
in the result is preserved. We conclude that, assuming a Milne- 
Table 3. Inversion's efficiency percentages in the weak field 
regime as function of the magnetic field strength. 



AB 


B <150G 


B <350G 


B <550G 


B <750G 


Noise 


lOOG 


75 % 


77 % 


81 % 


84 % 


No 


200G 


98 % 


99% 


99% 


99% 




lOOG 


60% 


57% 


61 % 


64 % 


5x10-" 


200G 


85 % 


87% 


89% 


90% 




lOOG 


56 % 


53 % 


54% 


57 % 


2x10-' 


200G 


76 % 


80% 


82% 


83 % 





Eddington model, the code developed in this work for the line of 
Mm at 553.7 nm, correctly inverts profiles in very different field 
strength regimes: from very weak magnetic fields (lower than 
100 G) up to 2 kG, with an eiTor bar of 200G with 89 % of prob- 
ability (300G with 96 % of probability) independent of the field 
strength. With such results we can assert that, even if the actual 
field strength cannot be measured from inversions with sufficient 
precision, at least the inversion algorithm can disentangle with- 
out any ambiguity the hG from the kG regimes of field strength. 



4. Observations 

We carried out two campaigns of observations at the T HEMIS 
telescope dLopez Ariste et al. 200"ot iGellv et al. 20061) in the 
multiline observation mode MTR. During the first campaign in 
June 2005, we searched for the spectral signature of the hyper- 
fine structure of the Mm atom among four candidate lines in 
the visible spectral region at wavelengths 539.4 nm, 547.0 nm, 
551.6 nm and 553.7 nm, respectively. From the data, we con- 
cluded that out of these four lines, only the one at 553.7 nm 
presented a detectable spectral signature of hps in the Stokes V 
profiles clear enough to be used for diagnostics. During the sec- 
ond campaign in June 2006 we focused on this line. In both 
campaigns the Fei doublet at 630.1 and 630.2 nm was ob- 
served simultaneously with the Mn i lines. Some relevant results 
concerning the importance of simultaneously observing spec- 
tral lines from both atomic species, as well as a brief descrip- 
tion of the data acquisitio n and reduction, has been presented in 
(iLopez Ariste et al' 20061) . 

In both campaigns we scanned the solar disk center and, in 
at least one day per campaign, we enjoyed excellent seeing con- 
ditions in the early morning. We have used the polarized signals 
levels of the Fe doublet (630.2 nm), observed simultaneously 
and co-spatially, to verify that around the scanned area there 
were no traces of active regions, nor enhanced magnetic areas: 
the Stokes V amplitudes of this magnetically sensitive line stay 
below 1% polarization at all times, indicative of mostly quiet 
photosphere, both network and intranetwork regions alike. The 
data sets consisted of temporal sequences of 100 modulation cy- 
cles in polarization. Each modulation cycle was made of 6 im- 
ages with exposure times of either 300 ms or 500 ms. The slit 
was kept fixed over the solar disk, up to residual seeing-induced 
image movements. In 2006 the new THEMIS correlation tracker 
was fully functional, thus further diminishing any residual move- 
ment. The total time of these sequences, around 8 minutes, is 
comparable to the lifetime of the granulation. The dimensions of 
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the slit spanned a region of 48 arc sec in two separate regions of 
15.5 arc sees each, with a masked region of 17 arc sec between 
them (this masked region is a requirement of the polarimetry 
mask used for the observation, which places the splitted, orthog- 
onal polarisation, beams 16 arcsec away one from the other in 
the direction along the slit). The slit width was kept at 0.5 arc 
sec resulting in a spectral resolution of 20 mA in agreement with 
the pixel spectral sampling of 10 mA. 

We reduced the data with the DeepStokes software tool cur- 
rently integrated to the THEMIS system facilities, modified to 
preserve spatial resolution. As delivered the DeepStokes code 
sums all images in the temporal series to increase the signal-to- 
noise ratio. Because of the long-time integration, image motion 
and blurring results in a degradation of the image quality. Trying 
to avoid such a degradation of the spatial resolution, otherwise 
attained in the individual images of the time series, we modified 
the code so that the addition of further images in the time series 
was not done regardless of a further degradation in the image 
quality. 

Specifically, in order to decide if the modulation cycle / + 1 
was to be added to the sum of the previous / modulation cycles, 
a cut along the sUt (at a pre-specified continuum wavelength of 
the intensity image) for the modulation cycle / H- 1 is compared 
to an analogous cut in the previous added / cycles. The number 
of points which differed in both such cuts in less than D was 
measured. When happened to be bigger than a pre-specified 
threshold A^o the ; + 1 cycle was deemed spatially coherent with 
the previous ones and its Stokes parameters added to the series, 
increasing the signal-to-noise ratio. If the / + 1 cycle failed to 
pass the test, it was assumed that the slit had moved and a new 
time series was started at point The sensitivity threshold D 
and the minimum number of cases A^o are provided by the user. 
We determined empirically their value from a posteriori consid- 
erations on attained signal-to-noise ratios and spatial resolutions, 
and the trade-off among both parameters. In the cases under anal- 
ysis in this work the achieved noise levels laid between 5x10"^ 
and 2x10"'* of the continuum level, for a spatial resolution better 
than 1 arc sec and approaching at times 0.6 arc secfl 

5. Results 

In Fig.|4]we show for illustration some of the observed Stokes I 
and V profiles (dotted lines) and the correspondent fits from the 
inversion (solid lines) overplotted. The panel covers the diversity 
of magnetic field strengths found in the quiet sun as seen through 
the Mn i lines: from the weak intensity regime (hG) to the strong 
regime (kG). 

5. 1 . The complete set of profiles 

Only 14% of the total amount of data had a S/N ratio high 
enough to be inverted reliably: around 500 profiles. In Fig.|5] we 
show the resulting magnetic field strength and longitudinal flux 
density distributions (defined as filling factor times the longitudi- 
nal field strength) from those inversions. For better description, 
we have fitted the histograms to an empirical distribution func- 
tion of the form N{x) = aexp(-x/b) + c, where a, b and c are 
positives constants and N(x) represents the number of cases. The 
introduction of a uniform distribution (the constant c) is required 



^ The spatial resolution was gauged by the crude and simple method 
of counting and measuring the size of brighter-than-average and darker- 
than-average regions as visible in the continuum cuts along the slit. 



to fit the flat part of the histograms, corresponding to intermedi- 
ate and high strength and flux density values. Otherwise, the left 
part of the histogram, associated to weak strength and flux den- 
sity values respectively, is well fitted with an exponential func- 
tion. In the case of the field strengh distribution, the fitting val- 
ues are a=672, b=51.3 G and c=8.6, while in the case of the flux 
density a=465, b=6.3 Mx/cm^ and c=1.7. The full widths at half 
maximum are respectively fwmhb = 54 G and fwmhh„ v - 4.6 
Mx/cm^, and the correspondent mean values are < B > - 635 G 
and < Flux > = 13 Mx/cm^. 

Although one expects m ostly maxwelhan-like probabihty 
distribution functions (PDF) (Dommguez Cerdefia et al. 2006bh 
for the field strength, because of the coarseness of the histogram 
bins we prefer to fit the data with an exponential for the weak 
fields plus a constant distribution for the strong fields. In spite 
of that empirical fit, the histogram shown in Fig. |5] should not 
be translated into that the most probable strength of the inverted 
profiles is zero G, but rather some value between < B < 75 G, 
included in the lower bin of the histogram. The main information 
the histograms show is that the weak field strength regime (hG) 
is more abundant or dominates the surface coverage in the quiet 
sun conditions (70%). Nevertheless, a considerable amount of 
cases, 30% of the inverted profiles, present magnetic strengths 
superior to 1000 G, i.e. in the strong field strength regime, and 
cannot be considered as absent. 

5.2. Network versus Intranetwork 

In order to study how the presence of network and intranetwork 
regions influence the results shown in the histograms of Fig.|5] 
we proceed to differentiate the profiles between these two re- 
gions. To this purpose, we have used the Fe i line (5538.52 nm), 
sitting next to the Mn line, to set a criterion for the classifi- 
cation of profiles into one or the other regions. We arbitrarily 
set the following condition for points to belong to the photo- 
spheric network: If the level of the normalized circular polarisa- 
tion V/I in the Fe i line exceeds a threshold value (TV) of the 
continuum level we consider the point as part of the network. 
All other points are considered intranetwork. We have initially 
used two threshold values: TV= 1x10"^^ and 2x10"^^ . In Fig.|6] 
we show both distributions (for field strength and longitudinal 
flux density) for each one of the regions. The solid line repre- 
sent a TV= 1x10"^ and the gray dotted line is for TV= 2x10"^. 
When using the upper threshold, TV- 2x10"^ corresponding to 
higher flux concentrations, the contribution of profiles in the 
weak strength regime in the network region has almost disap- 
peared while the number of profiles with strong field strengths 
has been just slightly modified. On the other hand, the use of the 
moderate threshold value 1x10"^ already identifies the profiles 
in the strong regime in the network region as well as a non neg- 
ligible contribution of fields with weak strengths. In the intranet- 
work set, with either one or the other threshold the general shape 
of the histogram remains unchanged. It is clear that a different 
threshold value able to radically change the general shapes of the 
histograms for both (network and intranetwork) regions would 
be unacceptably high and unrealistic. Therefore, and keeping in 
mind that small dependence of the network-assigned weak fields 
on the threshold value, we adopt in what follows the network 
and intranetwork definitions using the lower TV, with correspon- 
dent distributions function represented by the solid line in Fig. 
|6] To complete the justification on the criterion use for classifi- 
cation int o network and internetw ork, let us add that from previ- 
ous work dLopez Ariste et al. 2006 ) we know that the described 
classification, however arbitrary, successfully identifies most of 
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Fig. 4. Examples of inverted profiles. The dotted lines correspond to the observed data and the solid line to the respective fitted 
profile. The profiles were inverted without any smoothing, but in this plot a 2-pixel smoothing has been applied to observed data. 
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Fig. 5. Top panels: Histograms of the distribution of the magnetic field strength and the longitudinal magnetic flux represented in 
logarithmic scale. Bottom panels: The same histograms in normal scale. The dashed lines represents the fitted distribution functions 
N(x). The bin sizes are 75 G and 5 Mx/cm^ in field strength and flux, respectively. 



the network as identified visually over a magnetograph image. 
In terms of magnetic flux, the threshold value of 1x10"^, would 
correspond roughly to 10 Mx/cm-, and it is certainly biased to- 
wards including any flux concentration found in the intranetwork 
into the network, and conversely, any low flux region in the net- 
work will be included in the intranetwork dataset, a bias that we 
should remember in the conclusions. 

After classification, from a sample of 477 inverted cases, it 
was found that 70% belonged to the intranetwork region, while 
the resting 30% were from the network. As already noticed in 
the previous paragraph, the correspondent distributions for the 
magnetic field strength in the network region, represented by the 
solid line in the left panel in Fig.|6] clearly differs from the inter- 
network distribution and also from the complete dataset distri- 
bution: It presents an apparent growth in the number of cases in 
the strong field regime, in a striking bimodal distribution. On the 
other side, the intranetwork distribution keeps a similar shape to 
the one found before: The intranetwork appears to be predomi- 
nantly permeated by weak fields. 

The main conclusion arising from these separate distribu- 
tions is that the derived magnetic strength distribution from 
the inversions of the Mn i line in the visible range at 553 
nm, in the low-flux quiet sun or intranetwork region (tak- 
ing into account the known bias of our definition of in- 
tranetwork) is dominated in number of cases by weak field- 



strength fields. This is in accord with previous results ob- 
tained with the analysis by measurement in the infrared range 
(Khomenko et al. 2003: Lin 199l;lLin & Rimmele 19 99). in the 
visible (Keller et al. 1994.) and with the recently retrieved distri- 
butions fr om the very high spectral r esolution (0.3 arc sec) hin- 
ODE data (lOrozco Suarez et al. 20071) . Moreover, we found that 
the most probable expected strength is B < 100 G, close to the 
result obtained from simultaneous inversions from the infrared 
and from the visible recently published in Martinez Gonzalez et 
al. (2008). 

On the other hand, strong flux regions, mostly the photo- 
spheric network, rather show a bimodal distribution of fields 
with a contribution of weak field-strength fields comparable to 
that found in weak-flux regions, and an increased appea rance of 
kG fields jDonunguez Cerdena et al. 2006bl: ILites 2002h . 

In the table |4] we summarise these results in terms of inver- 
sion percentages for different magnetic strengths as seen in the 
histograms of the analysed data sets previously discussed. 

6. Conclusions 

In this work we contribute to explore the diagnostic capabili- 
ties of Mn lines with strong hyperfine coupling for the magnetic 
topology of the quiet sun. We built an inversion code for the Mn 
line at 553.7 nm formed in a Milne-Eddington atmosphere, and 
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Fig. 6. Histograms of network and intranetwork regions. The solid (dotted) line correspond to a threshold value of 1x10 (2x10 ^). 
In the network, the contribution of profiles with strength superior to 1 kG is considerably more important than in the intranetwork. 



Table 4. Comparison between the total set of inverted profiles 
and the Network-Intranetwork profiles. 





Total inverted 


Intranetwork 


Network 


Number of Profiles 


477 


339 


138 


Strong regime (kG) 


30% 


17% 


63% 


Weak regime (hG) 


70% 


83% 


37% 


<B> (G) 


635 


435 


1130 


< Flux > (Mx/cm^) 


12.8 


7.3 


26.5 



tested it under the conditions typical of the quiet sun in terms of 
signal amplitudes and noise levels. Such conditions are highly 
constraining: in most of the cases the magnetic field informa- 
tion is contained exclusively in the Stokes V profile, the linear 
polarization signals being under noise levels and the intensity 
profile being dominated by the non-magnetic portion of the res- 
olution element. Under such restrictive conditions the number of 
independent observables for the magnetic field is highly reduced 
jAsensio Ramos et al. 2007b and it sets a limit to the number of 
free parameters of our model. Previous works have shown that 
the Mn lines increased the number of observables thanks to the 
spectral features appearing under certain magnetic regimes. The 
inversion code we tested in this work takes advantage of it and al- 
lows us to improve over the simplistic weak-field approximation 
of previous works into a more sophisticated Milne-Eddington 
atmosphere and yet to conserve the possibility of disentangling 
strong from weak fields. 

Since most of the other parameters in the Milne-Eddington 
model, aside from field strength and flux, were not determined 
accurately under quiet-sun conditions, we doubt that any more 
sophisticated model can be used in the diagnostics at this point 
without the addition of further independent observables (e.g. 
with other spectral lines). However, more sophisticated models 
can shed Ught into implicit biases in such unrealistic picture of 
the solar photosphere as a M ilne-Eddington atmosphere is, e.g. 
(ISanchez Almeida et al. 20081) . 

We have payed special attention throughout this work to the 
numerical tests made to ascertain whether, under present obser- 
vational conditions of the quiet sun, the weak fields were disen- 
tangled from the strong fields for similar net fluxes. The answer 
is aflirmative and therefore we proceeded to apply the code to 
real data taken with the THEMIS telescope under good seeing 



conditions. Even in a telescope like THEMIS, focused on per- 
forming the most sensitive polarimetry possible, the acquisition 
of the data at the required signal-to-noise levels implied a trade- 
off in spatial resolution (seldomly better than 1 arc sec) that lim- 
its our results. We therefore constrained ourselves to separate the 
photospheric network from the intranetwork in our study, and 
restrained of identifying smaller scales, like granules and inter- 
granules. 

To identify network from intranetwork we applied a rather 
poor rule fixing at 10"^ the smallest amplitude of the Stokes V 
profile respect to the continuum intensity for a point to belong 
to the photospheric network. Such a rule obviously classifies as 
network any high flux concentration, mostly found over the net- 
work, but not totally absent from the intranetwork. Conversely 
it calls intranetwork any low flux concentrations, mostly found 
all over the intranetwork but not absent from the network either. 
With such bias in mind we decided to keep the names of network 
and intranetwork for the two classes because a small change in 
such threshold value did not change the general shape of the his- 
tograms, but also because of the following reasons. 

After inversion, the intranetwork points show a distribu- 
tion of field strengths whose dominant feature can be eas- 
ily fitted with an exponential tail. We use the noun tail be- 
cause an exponential distribution would result in a non-zero 
probability for null fields, something unacceptable for a vec- 
tor field. Such a distribution, showing that the quiet sun re- 
gions are mostly permeated of weak strength fields, would 
be in full agreement with previous results in both Han le and 
Zeeman effbct s (Manso Sainz et al. 2004; Lin & Rimmel e 19991: 
iKhomenko et al. 2003; .Martinez Gonzalez et al. 200 8) and re- 
veals a random or disorganised vector magnetic field at scales 
much smaller than present spatial resolution, if Hanle and 
Zeeman diagnostics are to be coherent with each other and with 
a field strength whose highes t probability remains below lOOG 
(ISanchez Almeida et al 2003h . 

The network points, on the other hand, show a bimodal dis- 
tribution: the exponential tail found in the intranetwork distri- 
bution is still present, peaking at field strengths below lOOG. 
But on top of that disorganised magnetic regime we observe the 
appearance of stronger fields, in the kiloGauss regime. This sec- 
ond distribution of fields would be in agreement with the pres- 
ence of one or several concentration mechanisms, able to sta- 
bilise, organise and perhaps amplify the statistical fluctuations 
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of the ubiquitous turbulent field in those places where photo- 
spheric dynamics or sheer accidents are able to maintain a tem- 
poral coheren ce bigger than the typical magnetic diffussivities 
jParker 1982b . One can heuristically think that the strong down 
draftin g plumes in the vertex of the co nvection cells are such 
places iNordlund et al. 1992HRast 20031) . thus defining the mag- 
netic photospheric network as a non-continuous succession of 
bright network points (Muller & Roudier 2003) where magnetic 
fields are organised in vertical structures and made stable by the 
downdrafting plasma. Other smal ler plumes found in mesogran- 
ular or granular scales ( Rast 2003 ) may also survive long enough 
to allow the statistical concentration of high fluxes, to be found 
here and there throughout the intranetwork and that, in our anal- 
ysis, would have been included in the network group. 

The cartoon just drawn above is in agreement with the his- 
tograms measured in the present work over the quiet sun. In 
particular, we are prone to the presence of an ubiquitous ran- 
dom field, as inferred from Hanle diagnostics, and whose sta- 
tistical fluctuations can be seen as Zeeman signatures over deep 
magnetograms, but otherwise highly disorganized at present res- 
olutions except when plasma movements create conditions, as 
around convection plumes, stable enough to concentrate and 
build up net magnetic flux. In any case, with the inversions per- 
formed here on the Mn i in the visible region of the spectra, 
at 553 nm, we have found that the distributions of the strength 
fields in the quiet Sun regions are dominated in occurrence by 
weak strength fields except in those areas related with the net- 
work. 
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